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Abstract. Motivated by a stochastic differential equation describing the dynamics of 
interfaces, we study the bifurcation behavior of a more general class of such equations. These 
equations are characterized by a 2-dimensional phase space (describing the position of the 
interface and an internal degree of freedom). The noise accounts for thermal fluctuations of 
such systems. 

The models considered show a saddle-node bifurcation and have furthermore homoclinic 
orbits, i.e., orbits leaving an unstable fixed point and returning to it. Such systems display 
intermittent behavior. The presence of noise combined with the topology of the phase space 
leads to unexpected behavior as a function of the bifurcation parameter, i.e. , of the driving force 
of the system. We explain this behavior using saddle point methods and considering global 
topological aspects of the problem. This then explains the non-monotonous force-velocity 
dependence of certain driven interfaces. 



PACS numbers: 05.45.-a, 05.10.Gg, 75.60.Ch 



Noise and Topology in Driven Systems - an Application to Interface Dynamics 



2 



1. Introduction 

In this paper, we consider a simplified, but general, description of driven dissipative systems 
described by two degrees of freedom in the presence of thermal noise. The theory applies 
to systems with two phases separated by a rigid moving domain wall (DW) with an internal 
degree of freedom, but it also describes general stochastic differential equations having a 
homoclinic saddle bifurcation. We will study in detail the behavior of such equations. 

The paper is written with two audiences in mind; those interested and familiar with 
dynamical systems in the presence of noise - and those more interested in physical 
applications. A short account has been given in [1 1]. 

1.1. Physical motivation 

A large variety of physical systems have interfaces separating different phases, with examples 
ranging from magnetic [12, 4, 23, 14] or ferroelectric [17, 18] domain walls, to growth 
surfaces [2, 10], contact lines [15]. The properties of an interface are well described at the 
macroscopic level by the competition between (i) the elasticity, which tends to minimize the 
interface length, and (ii) the local potential, whose valleys and hills deform the interface so as 
to minimize its total energy. 

The theory of disordered elastic systems [9, 8] allows one to determine their static and 
dynamical features {e.g., the roughness at equilibrium and the response to a field). Applying a 
force /, the interface can be driven to a non-equilibrium steady state. A crucial feature of the 
zero-temperature motion is the existence of a threshold force /^.^j^ below which the system is 
pinned. The system undergoes a depinning transition at / = f^^.-^^ and moves with a nonzero 
average velocity f for / > /^.^.j^. Close to the transition the velocity f ~ (/ — fa-itY 
characterized by a depinning exponent (3. In all these situations, the velocity is a monotonous 
function of the force (the more the interface is pulled, the faster it moves). Predictions of this 
theory are in very good agreement with experimental results, especially in the creep regime 
for interfaces in magnetic [12] or ferroelectric [18] films. 

In spite of this success, there are situations where the disordered elastic theory does 
not apply: for instance, one basic assumption is that the bulk properties of the system are 
summarized by the position of the interface alone. Here, we study the case where the position 
of the interface is coupled to an internal degree of freedom and we will show how this 
coupling affects the motion of the interface. An example is provided by domain walls in 
thin ferromagnetic films, where it is known that such an internal degree of freedom (a phase, 
to be detailed below) plays an important role. 

We reproduce in Fig. 1 experimental measurements of the mean velocity. It is puzzling 
that the velocity is not a monotone function of the force. 

The aim of this paper is to shed some light on this problem, by discussing a very 
simplified version of the system. We explain the general features of Figure 1 by two 
ingredients: first, by observing that there is a change of the topology of a typical evolution as 
a function of the driving force, and second, by taking into account temperature. While we will 
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Figure 1. The experimental velocity (dots) of an interface in a narrow nanowire, driven by 
a small external current, adapted from [16]. The Walker model (represented by the dashed 
line), which discards the pinning potential, does not reproduce the experimental results. The 
horizontal axis is the force (magnetic field (Oe)), the vertical axis is velocity (m/s). 



work with a simplified potential, we will gain some quite general insights on this and related 
problems. 

The experiments mentioned above come with physical models which describe the 
interaction between the phase and the position r of the wall. For the purposes of this 
paper, we will use the rigid wall approximation, [20, 13, 21, 22, 3, 5]: 

ad^r -dtip = f - cos(r) + r]^ , 

The external field / — cos(r) describes a constant "depinning" (or "tilt") force / and 
a "pinning" force — cos(r) deriving from a periodic potential. The damping coefficient 
a accounts for Gilbert dissipation. The effective thermal noise is a white noise with 
correlations [5] {'ni{t)Vjit')) = 2{hN)-^akBT5{t' - t) 5^^ where N = 2XA/a^ is the number 
of spins in the DW, where A is its cross-section, a the lattice spacing. Last, Kj_ is the 
anisotropy constant of the ferromagnetic medium. 

1.2. Mathematical motivations 

The study of (1) reveals that the system has a saddle-node bifurcation at / = /s = 1. 
Furthermore, for a range of fixed Kj_ one finds values of / for which the unstable manifold of 
the unstable fixed point is homoclinic.| 

I More precisely, one side of the unstable manifold is homoclinic, while the other goes to a second (stable) fixed 
point (H-^ and S in Figure 4 and Figure 5). 
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Saddle bifurcations have been discussed in many different contexts, and the influence 
of noise is well studied. Early papers are [6] and [1]. In those papers, the setup is that of 
intermittency in the presence of noise, with a discrete dynamical system of the form 

Xi+i = Xi-e + xj + ii + h{Xi) , (2) 

where G M, £ G M is the bifurcation parameter, and is some appropriate noise. The term 
h{x) describes a function which, e.g., vanishes in the neighborhood of x = but is such that 
orbits must eventually return to a neighborhood of 0. For this setting, under weak additional 
assumptions, one can study in quite some detail the invariant measure, and several salient 
features appear: 

• Orbits stay a very long time close to a; = and do fast excursions away from the 
neighborhood. 

• When the parameter e is positive, the deterministic system has a stable and an unstable 
fixed point (close to a; = -^y/e). The stochastic dynamics then helps the system to escape 
from the attracting fixed point (which is at ~ — \/?), but this may take a long time. 

In this paper, we discuss a similar scenario, but with some new features: We consider the 
parametrization f = 1 — e^. 

(i) There are two equations (and they are differential equations rather than iterations), with 
a saddle-node bifurcation at the value of the bifurcation parameter e = = 0. 

(ii) Close to £ = there is an ejj > for which the unstable manifold (of the unstable fixed 
point) returns to the unstable fixed point. § We will be interested first in what happens for 

(iii) We then discuss how the topological type of the orbits can change when the phase space 
is a torus. This will lead to a non-monotonous mean sojourn time near the unstable fixed 
point. 

The normal form of (1) is obtained by various rescalings, and a non-linear coordinate 
transformation. The deterministic part is given by 

dx = {ex + x'^) dt , 
dy ^ -ydt . 



While the deterministic part follows in a quite simple way, there is also a term appearing 
from the change-of-variables (the Ito term).|| This term takes the form —a^Qdt (in the x- 
component above, and a similar term for the y-component) with 

§ This is called a homoclinic connection. 

II We thank a referee for pointing out our oversight in not discussing this term in an earUer version of the paper. 
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when / = 1 — We will study (1) in the regime where £ > and o"^ <^ £. (The simulations 
of Figure 3 were done for e > 0.03 and cr^ < 8 ■ 10^^.) 

Therefore, we continue the discussion of the local equation near the fixed point with the 
more easily tractable form 

dx — {ex + x^) dt + cr^d^j. , 

dy ^-ydt + ayd^y , (3) 

and omit the Ito term. Here, d^^ and d^^^ describe the white noise, and the three parameters 
are e > and > 0, ay > 0. Adding a term h as in (2) on can achieve that globally 
the unstable manifold of x = y = returns to a; = y = for some small £h > when 
a — 0. We will tacitly assume that such a term is present. The phase space of (1) is the torus 
(r, (p) e [0, 27r) X [0, tt) and the unstable fixed point is at r = 0, </? = 0. We will argue in 
Sect. 5.2 that the term (Jyd^y can be omitted without changing the qualitative behavior of the 
problem. 

2. Results 

We first present the results from a physicist's perspective: 

In Figure 2 we illustrate the first two findings which appear because a second field cp 
comes into play: 

• The critical force, at which depinning initiates, moves from /g, ( = 1), to a lower value 
/h. Between /jj and /g the system is bistable: the velocity is either or strictly positive 
(see Figure 2). 

• The critical exponent of the velocity at depinning changes from | to "+oo ": the velocity 
grows like v l/\ log(/ - /h)|. 

The physical picture behind the bistable regime is the following: The position r 
represents the position of a particle in a tilted periodic potential. For / > /g this potential 
presents no local minima and the velocity is positive. For / < /g there are local minima that 
cannot be overcome in the absence of (p (this corresponds to the dashed curve of Figure 2). 
The phase p acts as an "energy store" for the position r. If r starts close to a local minimum, 
dissipation makes it end at the minimum and the steady velocity is zero (this is the lower 
branch of the bistable regime in Figure 2). On the other hand, if r starts far from a local 
minimum, the system reaches a stationary regime where (p helps r to cross the energy barriers 
between successive minima, by periodically borrowing and giving "kinetic energy" to r (this 
is the upper branch of Figure 2). 

Furthermore, if we introduce temperature, i.e., some extemal noise, then the force- 
velocity curve no longer presents any bistability (see Figure 3). This leads to a third 
observation: 

• The appearance of a third critical force, /-p. Note that for all (small) positive 
temperatures T > 0, the force-velocity curves actually cross at f = /t, as illustrated in 
Figure 3. 
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Figure 2. The two scenarios for the mean velocity: The green (dashed) curve shows v for the 
case of (4) where no internal degree of freedom is present, which corresponds to the motion 
of a particle on an inclined "washboard". The blue (solid) curve shows the velocity v for the 
case of (1) where the domain wall position r is coupled to tp. 




Figure 3. The velocity for (1) as a function of / and for several small values of the temperature 
T. Note that the curves cross at some value /-p. In the limit T J, the curves accumulate at 
/rp, while the deterministic equation {i.e., T = 0) leads to the blue curve. The parameters 
are e {0.03,0.1}, a £ {0.0001.0.0009}, a = 1/2 and = 6. The integrator was 
Euler-Maruyama, with a time step of 0.0003. We averaged over 512 samples. 



Our next result is a consequence of the periodicity of ther.h.s. of (1) in This periodicity 
is typical of domain walls in magnetic systems. The domain wall position is generically 
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coupled to an internal degree of freedom (for example a phase </?). ^ In Section 6 we will 
show how the periodicity influences the mean velocity 

• The mean velocity of the system (1) is a non-monotonous function with many maxima 
(depending on the values of a and K^). 

Remark. In early work [20], it was observed that, in the absence of pinning, i.e., because 
of the cosine in (1), increases up to a characteristic force (called the Walker force) 
above which the velocity decreases for a large range of / > see Figure 1. What we show 
is that the pinning potential leads to a very different scenario. 



3. A simple example 

Leaving for a moment (1) aside, we study first a much easier problem to familiarize the reader 
with our approach. We consider the problem of depinning from a periodic potential, but 
without the phase </?. The common underlying ingredients of such systems is the "pulling" of 
an interface by a force /. As the easiest example, we can consider the case of an "inclined 
washboard": 

dtr^ f - cos(r) , (4) 

where / is the constant force and r = r{t) e R is the position of the DW at time t. Clearly, 
the r.h.s. of (4) can vanish only if |/| < 1, and in that case every initial condition Tq = r(0) 
will, as time evolves, converge to one of the values = arccos(/) + 27rn, with n any integer. 
In this case, we say that the potential is pinning. On the other hand, when |/| > 1, there is no 
fixed point for (4) and r{t) will increase or decrease indefinitely. In fact, one can check that, 
for / > 1 and r(0) = the solution of (4) is: 

/tan (1^77^) V<F^ 
r{t) — 2 arctan I f + 1 

whose derivative is a periodic function of t with period p — 27r/-\//^ — 1. Therefore, 
r{p) — r(0) = 27r, and the limit is 

lim!:W=^ = 0(V7^). 

t^oo t 27r J J 

In other words, the mean displacement, which we call the velocity v, is given by ?;(/) = 
a//^ — l/(27r). Thus, for the simple case of (4) the well-known result is that near the 
depinning transition, the velocity grows like >//g — / where /s = 1 in our simple example. 
Furthermore, the velocity is obviously a monotone function of /: The harder one pulls, the 
faster one advances. 

f Although this coupling is well known in the magnetic DW community [13], it has to our knowledge always 
been discarded in interface physics. 





Figure 4. Phase portrait of the evolution (1) in the (r, 93) plane, for £g = < e < Ejj. 
Boundaries between different regions are in red. The attracting limit cycle is shown as a 
yellow line, while the repulsive limit cycle is shown as a yellow dashed line. In green, the 
basin of attraction of the stable fixed point S ; in blue, the basin of attraction of the stable 
limit cycle, either from below (light blue) or from above (dark blue). 



4. The coordinates of the problem at zero temperature 

We now study the special case of (1) when the noise terms are absent 

adtr -dt<^ = f - V'{r) , 
adt^ + dtr = - sin(2y?) , 

where V'{r) = cos(r). 

When Kj_ is very large, Lp will be very close to mod vr, and then the system reduces 
to the washboard model (4). However, for smaller Kj_, the phase matters and this is the 
case we want to study now. A redefinition of / = 1 — brings the problem (5) to the more 
convenient form 

adtr - dtif = -e^ + (1 - cos(r)) , 

(6) 

ad^if + d^r = —^Kj_ sm{2ip) . 

The phase space of this equation is the torus {r,Lp) G [0, 27r) x [0,7r). For the following 
discussion the reader is referred to Figure 4 where the torus is drawn in the plane with the 
horizontal axis corresponding to r and the vertical corresponding to (/?. A three-dimensional 
rendering is shown in Figure 5. 

It is easily verified that (6) is invariant under the symmetry: r — )■ —r, cp — )■ —tp + ti/2, 
t — >■ — t. This makes the phase space centrally symmetric, but we will not make use of this 
property in the analysis. 
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Figure 5. The unstable manifold (purple) of the (yellow) hyperbolic fixed point winds 
around the torus once (counterclockwise) and ends at the fixed point . In green (behind) the 
same for the other fixed point iJj- The stable fixed point is at the end of the blue "tail", and 
the unstable at the end of the orange tail. 

We will consider only values of e'^ > 0, K^^ > 0. For the simulations, we took a = ^. 
Under these assumptions, the local structure of this equation is characterized by 4 fixed points 
of the form (0, ir^) and {tt/2, ir^), where 

= arccos(l — e^) . 

The stability of the 4 fixed points is as follows (for e > 0): 

• Hi = {0,r^) and H2 = (vr/2, — r^) are hyperbolic (with one stable and one unstable 
direction), 

• S = {0, — r^) is stable, 

• U = (11/2, r^) is unstable. 

For e = we have = and the corresponding pairs of fixed points collide, leading to a 
single fixed point with one direction stable, and the other stable-unstable. Thus, at e = the 
fixed points S and (resp. U and H2) collide; we are in the presence of a typical saddle-node 
bifurcation. 

5. General discussion for the case of non-zero temperature 

Apart from its interest as a physics problem, the equations under study are a nice example 
of the interplay of homoclinic orbits, collision of a saddle-node, and the influence of noise. 
While any combination of two of the three phenomena is amply discussed in the literature. 
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[19], as far as we know, the combination of all three seems to be new. In particular, as we 
shall show, the system will have a "phase transition" as the noise goes to zero, which occurs 
neither at the homoclinic point, nor at the collision of the saddle-node, but at a well-defined 
intermediate point. The present section will derive this in a general form. 

5.1. The one-variable case 

In very early work, Risken [19], considered the problem 

9fr = —^d^r — £ + (1 — cos(r)) . 
If we write it as a first order system, we have 

d^v = —^v — £ + (1 — cos(x)) . 

The phase space for this system is {x,v) E [0,2n) x M. There are now only two fixed 
points: v = 0, x = = ±arccos(l — e). So, the system is really quite different from our 
model. However, two of its main features remain and they can be discussed in the spirit of 
(1): Locally, there are two fixed points: One is stable {v, x) = (0, x^) and the other (0, —x^) 
is hyperbolic. Again, for £ = there is collision of the two fixed points (a saddle-node 
bifurcation). On the other hand, there is a value £h of ^ (^lot the same as in our model) 
depending on 7 for which we have a homoclinic connection). 

5.2. The 2-variable case 

We consider again (1), but change coordinates immediately to a normal form. Furthermore, 
for the purpose of the discussion in this section, it is irrelevant that the natural phase space is 
the torus. In fact, it suffices to consider a local coordinate system near the saddle-node. The 
global aspects only have to do with the "reinjection" [6]. 

In a local coordinate system where the hyperbolic fixed point H^^ is at the origin, up to 
terms of higher order, and neglecting the Ito term, as discussed in Section 1.2, the system can 
be written in the form 

dx = [ex + x^) dt + crd^ , 
dy = -y + (72d^2 ■ 

Here, d^ and d^2 describe white noise, and the three parameters are £ > 0, cr > 0, and (T2 > 0. 

We will omit the noise term in the (stable) y direction because it would only induce 
a fluctuation in the "arrival time", but has no influence on the escape rate to the basin of 
attraction of the fixed point (— £, 0). However, the noise in x direction is essential for our 
discussion. 

There is one more, crucial, assumption: For some £h > (when a = 0) the unstable 
manifold of the fixed point {x, y) = (0, 0) (in the positive direction) is homoclinic, that is, it 
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returnsto (0,0). Furthermore, for e < the unstable manifold is moved to the right (positive 
x). See Figure 6. We also assume that this unstable manifold is transversally stable, that is, 
nearby orbits are attracted to it, as illustrated in Figure 6. Such behavior can be obtained if 
in (7) we add some nonlinear terms which bend the unstable manifold of (0, 0) as shown in 
Figure 6. We assume in the following that (7) has been modified accordingly, without changing 
the vector field near (0, 0). 

Proposition 5.1. Under the above assumptions, there is a constant A> Q such that the mean 
velocity of the system ( 7) has a phase transition at a point e-p and, for small e and large e^/ a^, 
this transition happens at £rp close to the solution of 

e-6{A- {eje - l)f = 0. 

(This solution lies in the interval (eg = 0, e-^).) 

The essential thing here is that £fj > 0. When a > the following happens: If we start 
at some point of the unstable manifold, and evolve with the noisy evolution, the orbits come 
back, for e between £3 = ^H' a basically Gaussian distribution around the unstable 
manifold, see Figure 7. 

At this point, an intriguing competition between two phenomena occurs. On one hand, 
because a > 0, some orbits (those on the "inside" of the homoclinic loop in Figure 7) are 
accelerated by the noise, since they avoid the close passage by the fixed point (0, 0). On the 
other hand, those which return to (0, 0) on the side of a; < fall into the basin of attraction 
of the fixed point (— e, 0) and they will need a long time to escape from that basin. This 
phenomenon has been studied long ago under the name of "intermittency in the presence of 
noise" [6]. In that case, it was always (rightly) assumed that the reinjection density is close 
to uniform across the basin. In the case at hand, the novel problem is that the probability to 
fall into the basin of attraction of the point (— e, 0) decreases as e decreases from e^. This 
is because the center of the probability distribution of orbits moves away from the basin as e 
decreases, see Figure 8. 

To quantify this phenomenon, we assume that to lowest order, the unstable manifold is 
moved by an amount A • (sjj — e) in the positive direction, i.e., ^4 > 0. The potential along 
the X-axis is shown in Figure 9. 

We next ask, for a fixed e E (0, ejj) ^ fixed x E [0, e] how long the stochastic process 
(7), starting at x, needs to escape to the right (to +00). We will neglect the y coordinate in 
this estimate. As is well known, and for example done in detail in Section 3 of [6], see also 
[7], this time is given by the Green's function of the differential operator G, (7): 

with Dirichlet boundary condition at x — +00. The expected time r{x) to escape from x is 
then given by 



r[x)^— dze-'^(^) / dwe'*^^), (9) 

(^^ Jx J -00 
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Figure 6. The phase portrait at e = £ 




X 



Figure 7. The phase portrait at £ < £jj. The unstable manifold of (0, 0) has moved to the right 
by a distance D A - {s^ — s) . 



with the potential h = —V given by: 



2 fez^ 



cr2 



The integral (9) can be estimated as in [6]. First one changes variables io u — z -\- w and 
V — z — w and finds 



2 



r{-^)^- j du dvexpi-^{-v' + {— + —)v)). (10) 
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X 



Figure 8. The same phase portrait as in Figure 7. Superposed is the (Gaussian) distribution of 
noisy orbits returning along the unstable manifold of (0, 0). Note that the relation between D 
and the width a of the distribution detemiines how frequently a noisy orbit will fall onto the 
stable (left) side of the y axis. 




Figure 9. The typical shape of the effective local potential V{x) = —ex'^/2 — near 
a; = 0. Note that the depth of the potential and its width depend on e. 



(Pushing the integration Umit to x — — oo is justified by the fact that anyway, most of the time 
is spent near x — —e.) The u integration can be done expUcitly and leads to 

We rescale hy v = ew and thus find 

-(-~) = — I du>(-)V^exp (-^(- --)). (12) 
Using the saddle-point approximation (the critical point is at w = 1) this integral behaves, for 



Noise and Topology in Driven Systems - an Application to Interface Dynamics 



14 



large e^/a'^, and neglecting the prefactor in front of the exponential, as. 



t(— oo) ~ exp 




) 



On the other hand, as illustrated in Figure 8, the probability to reach a point x < is 
proportional to exp(— const.(|a;| + F^/a'^), where the constant F depends on certain global 
aspects of the problem, such as the length of the (almost) homoclinic loop. This just estimates 
how much probability leaks to the "wrong", i.e., left side of the unstable manifold of (0, 0). 
We will continue the discussion by assuming all the constants to be 1. A rescaling of the 
variables would eliminate an arbitrary constant anyway. In particular, the average time to 
leave the trap (say, between x — —e and x = 0) is then given approximately by 



For fixed e^, this polynomial has exactly one real root — £t(^) ^h) which lies in (0, £h)- 
This is the point where the behavior will switch over It is the point which corresponds to /p 
in Figure 3. 

6. Global topological aspects 

After having neglected the torus structure of the problem, we reinstate it in the current section. 
If we want to perform a global study of the system in the parameters e and Kj^ we have to take 
into account that the phase space of (6) (or (1)) is a torus. Thus there can (and do) exist several 
topologically different ways in which a homoclinic orbit can form. They can be indexed by 
two (non-negative) integers g and 'ip which count how many turns of 2n the variable r resp. 2(p 
will undergo as one moves from the fixed point to reach it again through the homoclinic 
loop, and denote by W(f?, ^p) the index of the orbit. 

In the space of e and Kj_, the picture which emerges numerically is shown in Figure 10. 
For each of the curves in Figure 10 we show one example in Figure 11. 

It is now easy to explain the non-monotonicity of the mean velocity for (1), as illustrated 
in Figure 12. Fixing a Kj^ (say Kj^ — 3.5 in Figure 10) and varying the pulling force f — 1 — e 
for e from to 0.3 we first cross the W(l, 1) curve and then the W(l, 0) curve. This leads to 
Figure 12. Note that, in accordance with the theory of Section 5, as a function of the noise 
(temperature), both bumps are filled with details which look as in Figure 3. In particular, there 
will be a special bifurcation point of the form frj. for both of them. 

In Figure 12, we illustrate the cases W(l, 0) and yV{l, 1). The reader should notice that 
for every pair (g^^jj) which is realized, for fixed Kj_, there will be a window W^,^ of values 




(13) 



Here, we have used that, to lowest order, D = A ■ {e-^ — e). 
Consider now the polynomial in (13). It can be written as 




{e-6{A.{e^/e-l)y). 
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Figure 10. The locus of some homoclinic connections in the e, plane. Shown are 3 such 
curves with winding number in r equal to 1, and in equal to 0, 1, 2. There are infinitely many 
such curves. 




Figure 11. From left to right: 3 homoclinic orbits, of type yV(l, 0), W(l, 1), and 2), 
respectively. 

of / around fuiQ,'^) which is like the case we discussed in detail. Mutatis mutandis our 
analysis will apply immediately to all these cases, as soon as the general hypotheses (about 
the transverse stability of the homoclinic orbit and the motion of the return as a function of /) 
are satisfied. 

The physical interpretation of the non-monotonous behavior of the velocity is the 
following. We have seen previously (see Section 2) that "helps" r to cross the barriers 
between the local minima of the potential where it lives. Doing so, ip oscillates in its own 
local minimum (this corresponds to the first bump in Figure 12). However, for larger values 
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Figure 12. A schematic illustration of the velocity as a function of force. First the topological 
type yV(l, 0) leads to a monotone increase of the velocity. But at some point, the winding 
number for changes, and the velocity drops to 0. Several other winding numbers, not shown, 
could occur before the topological case yV(l, 1) sets in. 



of the force, if will itself cross the barriers of its potential and dissipate so much energy that it 
cannot help r anymore (in the phase space picture of Figure 4, this corresponds to a collision 
between the attractive and repulsive limit cycle). There is a whole regime of force where no 
limit cycle exists (this is the flat region between the two bumps in Figure 12). It is only for 
larger values of / that a stationary regime appears where both Lp and r can cooperate and 
display non-zero mean velocity (this corresponds to the second bump in Figure 12). 
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